tridiagonal systems in matlab